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Magnetization switching due to a current-pulse in symmetric and asymmetric spin valves is studied 
theoretically within the macrospin model. The switching process and the corresponding switching 
parameters are shown to depend significantly on the pulse duration and also on the interplay of the 
torques due to spin transfer and external magnetic field. This interplay leads to peculiar features in 
the corresponding phase diagram. These features in standard spin valves, where the spin transfer 
torque stabilizes one of the magnetic configurations (either parallel or antiparallel) and destabilizes 
the opposite one, differ from those in nonstandard (asymmetric) spin valves, where both collinear 
configurations are stable for one current orientation and unstable for the opposite one. Following 
this we propose a scheme of ultrafast current-induced switching in nonstandard spin valves, based 
on a sequence of two current pulses. 

PACS numbers: 67.30.hj,75.60. Jk,75.70.Cn 



I. INTRODUCTION 

The possibility of current-induced magnetic switching 
(CIMS) in multilayer structures has been introduced by 
Slonczewsk^i and Bcrgci— . They pointed out that spin 
polarized current passing through a multilayer consist- 
ing of two ferromagnetic layers separated by a nonmag- 
netic spacer can exert a torque on the local magnetic mo- 
ments in ferromagnetic layers. In case of commonly used 
(standard) spin valves, this torque can switch the system 
between parallel (P) and antiparallel (AP) alignments 
of the layers' magnetic moments. Indeed, the CIMS 
has been confirmed experimentally in many spin valve 
structures^iii^i^. 

In order to meet the technological demands for ap- 
plications, designing of bistable magnetic spin valve de- 
vices that could be fast switched by electric current only 
is highly desired. Accordingly, several subtle switching 
schemes aimed at speeding up the switching process and 
lowering the energy costs have been developed for vari- 
ous magnetic multilayer structures. It has been shown 
that proper optimization of the current pulse parameters 
(amplitude, duration, frequency, etc.) can finally result 
in ultrafast, single-step, switchingii^. 

Current-induced dynamics in spin valves is governed 
by spin-transfer torque (STT), particularly by its mag- 
nitude and dependence on the angle between magnetiza- 
tion vectors of the reference and sensing magnetic layers. 
Such an angular variation depends mainly on the type 
of electron transport through the spin valve. Recently, 
STT in the diffusive transport regime draws more atten- 
tion. It has been demonstrated^, that spin diffusion in 
various types of polarizing layer can strongly influence 
the current-induced behavior of spin valves. To calculate 
STT in the diffusive transport limit, one can consider the 
modeUS. that extends the Valet-Fert descriptionii to ar- 
bitrary magnetic configuration and allows to study STT 
as a function of layers' thicknesses, type of materials, spin 



asymmetries, etc, which can be directly compared to ex- 
perimental observations^^. In the case of diffusive trans- 
port, spin-diffusion length has to be much longer than the 
mean free path. However, it has been showni^ that dif- 
fusive approach is valid also for spin diffusion lengths of 
the order of mean free paths. Moreover, the approach de- 
scribes relatively well experimental results even for layer 
thicknesses comparable to the corresponding mean free 
paths. 

One of the most pronounced manifestation of the 
importance of spin-diffusion length has been shown 
for asymmetric spin valvesiS, where the non-standard 
(wavy-like) angular dependence of STT leads to the 
current-induced precessional regimes in zero magnetic 
fieldi^ii^ii^. This prediction has been later confirmed 
experimentall y^^! . The wavy-like STT vanishes not 
only in the collinear configurations - P and AP ones - but 
also in a certain non-coUinear configuration. Depending 
on the current direction, the collinear configurations are 
either both stable or both unstable. The first case is of 
importance for the current-induced zero-field microwave 
excitations, while the later one would be desired for fur- 
ther stabilization of the collinear magnetic configurations 
in memory devices against thermal or current fluctua- 
tions. Since the problem of switching in non-standard 
spin valves has not been addressed yet, here we present 
a comprehensive study of the CIMS in such systems due 
to a current pulse of finite duration. In the following, the 
spin valves with wavy-like STT will be rcfcreed to inter- 
changeably as nonstandard or asymmetric ones, whereas 
the spin valves with STT vanishing only in collinear con- 
figurations will be referred to as standard or symmetric 
ones. 

In this paper, we investigate current-pulse-induced 
magnetization dynamics in three-layer spin valves 
X/Cu(10)/Py(8) sandwiched between semi-infinite Cu 
leads. The numbers in brackets correspond to the lay- 
ers' thicknesses in nanometers. The Py(8) stands for the 
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permalloy sensing layer, which magnetization direction 
can freely rotate upon applied magnetic field and/or spin 
polarized current, whereas X denotes the reference layer, 
which magnetization is considered to be fixed and not 
influenced by external magnetic fields and electric cur- 
rent. Assuming these conditions, we have considered two 
different types of the reference layer: X = Py(20), and 
X = Co(8). Taking into account the above introduced 
notation, the former spin valves, Py/Cu/Py, we will re- 
ferred to as standard or symmetric, whereas the latter 
spin valves, Co/Cu/Py, we will be referred to as non- 
standard or asymmetric ones. 

The main motivation of this paper is to provide a sys- 
tematic study of the CIMS, based on the macrospin sim- 
ulations. We consider a rectangular current pulse char- 
acterized by the current amplitude and pulse duration. 
We have found that a proper choice of the pulse param- 
eters leads to fast switching with relatively low energy 
costs. In the case of non-standard spin valves, the double- 
pulse scheme, which allows reliable switching between the 
coUinear magnetic configurations and significant shorten- 
ing of the overall switching time, is proposed. 

The paper is organized as follows. In section 2 we de- 
scribe the macrospin model. Numerical analysis for both 
standard and non-standard (asymmetric) spin valve is 
described in section 3, where also a double-pulse switch- 
ing scheme for an asymmetric spin valve is described. 
Section 4 includes summary and final conclusions. 



II. MACROSPIN MODEL 

Time resolved imaging of CIMS showed that inhomo- 
geneous spatio-temporal magnetization evolution takes 
place during the reversali^. Thus, more sophisticated 
models, involving detailed micromagnctic description, 
have to be adopted for a better qualitative description22,. 
However, the macrospin model - particularly for spin 
valves in nano-size range ~ provides a sufficient quantita- 
tive description of the current-induced dynamics, which 
is in agreement with many experiments on standard spin 
valves^. A relatively good qualitative agreement has 
been reached also for non-standard (asymmetric) spin 
valves, where some predictions based on the macrospin 
modeli^ii^'^ have been confirmed experimentall y^^i^^ . 

Time evolution of the unit vector s = {sx, Sy,Sz) along 
the net spin moment of the sensing layer is described by 
the generalized Landau-Lifshitz-Gilbert equation. 
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where 7g is the gyromagnetic ratio, /iq is the magnetic 
vacuum permeability, Ms is the saturation magnetization 
and d is the sensing layer thickness. The Gilbert damping 
parameter a is assumed to be constant, a = 0.01. The 
effective field Hcs includes contributions from the exter- 
nal magnetic field (i?cxt)i uniaxial magnetic anisotropy 
(^fani), demagnetization field (JJdom) and the thermal 
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FIG. 1: Angular dependence of the spin transfer torque 
components in units of hi /\e\ acting on th Py(8) in the 
Cu/X/Cu(10)/Py(8)/Cu spin valves for the standard case 
with X = Py(20) (dashed lines) and nonstandard case (wavy- 
like torque) with X — Co(8) (solid lines), (a) In-plane com- 
ponent Tg, and (b) out-of-plane component r^. The other 
parameters as in Refs [lollla ]. 



field (i?th); Hcs = -Hcxtez - i?ani {s ■ e^) + Hdcm + 
Hth, where is the unit vector along the axis z which is 
parallel to the in-plane magnetic easy axis. By definition, 
the external field is positive when it is pointing along the 
negative z axis. The demagnetization field corresponds 
to the sensing layer of an elliptical shape with the major 
and minor axis of 130 and 60 nanometers, respectively, 
and thickness of 8 nanometers. The magnetic easy axis 
is assumed to be along the longer axis of the ellipse, and 
Hani = 100.5 Oe. 

The thermal field Hth contributing to Hcs is a stochas- 
tic field of statistical properties: {Hth,i{t)) = and 
(Hth,,(t)i^th,j(i')> = 2VS,jS{t~t'), where i, j e {x,y,z} 
The strength of the thermal fluctuations is given by the 
parameter V = 2kBT/[Ms{l + a^)] derived according to 
the fluctuation-dissipation relation2ii2^, with T being the 
temperature. 

As concerns the STT, we consider both the in-plane 
and out-of-plane components, t = Te + t^, where Tg = 
al s X {s X S) and = bl s x S. The vector S points 
along the net spin of the reference layer X, S = Bz, and 
is assumed to be constant in time (current does not ex- 
cite magnetic moment of the reference layer). Current 
density / is defined as positive when current flows from 
the reference layer towards the sensing one. Finally, the 
angular dependence of the parameters a and b has been 
calculated in the diffusive transport limit^^. The result- 
ing in-plane and out-of-plane torques for the spin valves 
X/Cu(10)/Py(8) with X = Py(20) and Co(8) are shown 
in Fig. 1 as a function of the angle 9 (azimuthal angle) be- 
tween the magnetic moments of the reference and sensing 
layers {s ■ S = cos 9). 



III. NUMERICAL RESULTS 

Here we study switching in the X/Cu(10)/Py(8) valves 
from the P to AP state due to a current pulse. Numerical 
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solution of the LLG equation has been performed using 
the Heun Scheme»2^ with an auto-adaptive time step. For 
the initial configuration we assume a biased state with 
9{t = 0) = 1° and cp{t = 0) = tt/2, where (p is the po- 
lar angle describing orientation of the vector s from the x 
axis parallel to the current flow. The current pulse i{t) of 
constant current density / and duration tp is applied at 
t = 0, i{t) = I[eit)~eit-tp)], where e(x) = 1 for a; > 
and 9(a;) = for x < 0. A successful switching event 
with the corresponding switching time tg is counted when 
s7(ts) < —0.99, where s7(i) is the exponentially weighted 
moving average^i, s7(t) = 77 Szit) -|- (1 - ??)sl(i - At), 
At is the integration step, and the weighting parameter 
?/ = 0.1. The moving average s7 is calculated for time 
t > t' when Sz{t') reaches the value of —0.9; otherwise 
s7(t) ~ Sz{t). From the experimental point of view it is 
more convenient to reformulate the switching condition in 
terms of the magnetoresistance. We note that this holds 
only for the standard spin valves, in which the magne- 
toresistance is a monotonic function of the 9 angl o^^''^^ . 
In asymmetric spin valves the magnetoresistance can be 
a non-monotonic function of 6^^^, and therefore direct 
calculation of the magnetoresistance is then needed. 



A. Dynamics in a Py/Cu/Py spin valve 

Let us study first the standard spin valve, X=Py(20). 
The in-plane and out-of-plane components of the STT 
acting on the Py(8) sensing layer show sine-like angular 
dependences, see dashed lines in Fig. 1. Switching time 
as a function of the pulse duration tp and reduced cur- 
rent density I/Iq {Iq = 10*Acm~^) is shown in Fig. 2(a). 
Here, we consider zero temperature limit and the exter- 
nal magnetic field is set to zero. Two different regions 
in the switching diagram can be distinguished. First, 
the white non-switching region is observed for short cur- 
rent pulses and low current densities. In this region, the 
energy gain due to STT does not overcome the Gilbert 
damping and system stays in the initial local magnetic 
energy minimum. The second region corresponds to suc- 
cessful switching to the AP state. The switching time 
ts is shown in the color scale. The ts decreases non- 
monotonously with increasing current density. The most 
bright area corresponds to the ultra-fast switching, in 
which the spin reaches the AP configuration before the 
current pulse ends {ts < tp). In such a case the switching 
is realized in a single ultra- fast step after a half precession 
around the x axis; see Fig. 2(b) right. 

The boundary between non-switching and ultra-fast 
switching regions develops into a ripple structure. In 
this region, the energy gain due to spin-transfer leads to 
a retarded switching, where the switching time tg > tp, 
see Fig. 2(b) left. The switching for t > is accompanied 
with a ringing, where the spin relaxes to the AP state due 
to energy dissipation via the Gilbert damping only. Such 
dissipation, however, is rather slow and therefore the re- 
tarded switching is much more slower than the ultra-fast 
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FIG. 2: (Color online) Current-pulse driven switching in the 
Py(20)/Cu(10)/Py(8) spin valve in the absence of magnetic 
field, (a) Switching time as a function of the pulse duration tp 
and reduced current density I/Io, lo = 10*Acm~^. (b) Tem- 
poral evolution of the spin component under rectangular 
current pulse (dashed lines) of amplitude I /lo = 1.4 and du- 
ration tp = 1.2 ns {retarded switching, left) and = 2 ns {fast 
switching, right), corresponding to the points marked in the 
switching diagram (a). 



To speed up switching from P to AP state one may 
consider a negative external magnetic field. On the other 
hand, a positive magnetic field exceeding the anisotropy 
field leads to commonly observed steady-state out-of- 
plane precessional (OPP) modes^i^i^, which are the re- 
sult of the energy balance between Gilbert damping and 
the energy gain due to spin-transfer. 

Our analysis shows that for positive external magnetic 
field, the continuous switching region in the diagram 
shown in Fig. 2(a) splits into an non-compact current 
dependent stripe structure. In Fig. 3(a) we show the 
switching diagram for iJoxt = 200 Oe. The switching re- 
gions alternate with the stripes where the spin transfer 
induces the OPP regime. Since, the current pulse is fi- 
nite, the final state depends on the actual spin state at 
the t = tp, which falls into the basin of attraction either 
of P or AP state. This is shown in Fig. 3(b), where two 
switching events under the current pulses of the same 
amplitude (J = 2.75 Jq) and different pulse duration are 
driven via the OPP regime. We note that in the switch- 
ing regions the spin dynamics is similar to the zero-field 
switching discussed above. 
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To elucidate the stripe structure, we plotted in 
Fig. 3(c) map of the final spin states as a function of the 
initial spin position Sq = s{t = 0), assuming constant 
current amplitude and pulse duration tp — > cx). The gray 
(black) regions correspond to the initial spin position, 
which results in the final OPP (AP) regime. Comparing 
the maps calculated for two different current densities 
(/ = 2.75/0 and / = 3.00/0), one concludes that the 
dynamical phase portrait depends rather strongly on the 
current density. In other words, current-driven dynamics 
from the same initial state can develop to different final 
states. Further increase of external magnetic field leads 
to shrinking of the P— >AP switching stripes. For fields 
much larger than the coercive field the switching stripes 
disappear so the OPP regime remains only. We note 
that the initial spin position (close to the P configura- 
tion) assumed in this paper (except Fig. 3(c)) is denoted 
in Fig. 3(c) by the circles. 

The sharp stripe structure is a result of deterministic 
dynamics and fixed initial condition. When a distribu- 
tion of initial configurations is taken into account, some 
smearing of the border between the stripes is observed 
(not shown) . The boundaries are also smeared when non- 
zero temperature is considered. In Fig. 3(d) we show 
the switching probabilities as a function of the current 
pulse density, calculated for pulse duration tp = 3 ns at 
T = 4.2K and T = 77K, and for fixed initial configura- 
tion. The statistics has been calculated from 10^ events 
for each value of the current density. The switching prob- 
ability follows the stripe structure in the zero tempera- 
ture limit, and decreases with decreasing current ampli- 
tude. For T = 77K, the probability is lowered by the 
factor of about 3 and the peaks broaden. From this fol- 
lows that smearing of the boundaries between switching 
and precessional regions increases with increasing tem- 
perature. In addition, positions of the peaks are shifted, 
which reflects the fact that thermal fluctuations act on 
the spin like an additional torque, and non-linearly influ- 
ence the spin dynamics. 



B. Dynamics in a Co/Cu/Py spin valve 

The Co(8)/Cu(10)/Py(8) spin valve exhibits non- 
standard STT acting on the Py(8) layer. Due to the 
wavy-like dependence of the STT, shown by the solid 
lines in Fig. 1, positive current stabilizes both the P and 
AP configurations. A negative current, in turn, desta- 
bilizes both the coUincar configurations. This charac- 
teristic property of the wavy torque raises the question, 
whether it is possible to switch an asymmetric spin valve 
between P and AP states without the need of an external 
magnetic field. 

In Fig. 4(a) we show the switching diagram from P 
to AP state under a rectangular current pulse. Here 
one may distinguish four characteristic switching regions. 
First region, denoted by (i), corresponds to low current 
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FIG. 3: (Color online) Effect of the external magnetic field 
/^cxt = 200 Oe on the current-pulse driven dynamics in the 
Py(20)/Cu(10)/Py(8) spin valve, (a) Switching time as a 
function of reduced current density ///o (/o = lO^Acm"^) 
and pulse duration tp. (b) Temporal evolution of the spin 
component under current pulses (dashed lines) marked in the 
switching diagram (a), corresponding to the amplitude / = 
2.75 If) and durations tp — 2.56 ns (upper part) and tp — 2.7 ns 
(lower part), (c) The maps of the final states as a function 
of initial spin bias for current pulse of / = 2.75 /o (upper 
part) and I — 3.0 /o (lower part), and for tp 00. The spin 
dynamics initialized from the points inside the gray (black) 
areas leads finally to OPP regime (AP state). Circles denote 
the initial spin bias used in the simulations, (d) Thermally 
assisted switching probability. Paw, from P to AP state, driven 
by a 3 ns rectangular current pulse, calculated as a function 
of current density at T = 4.2 K (solid line) and T = 77K 
(dotted line). 



amplitudes and/or short pulses, where switching does not 
take place. The non-compact region (ii) of rib-like struc- 
ture borders the non-switching region and comprises rela- 
tively short pulses leading to the fast switching processes. 
In region (iii) the P/AP bistability of the final states is 
observed. Finally, in the region (iv) the final state of the 
dynamics depends on the applied current density, result- 
ing in the band-like structure. The structure contains 
regions with final P state, which regularly alternate with 
the regions of final AP state. 

In order to explain the complex diagram structure. 
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let us study currcnt-pulse-induced dynamics due to the 
6 ns pulse of amplitude / = — 3.45/o at zero temper- 
ature. The temporal dependence of the spin compo- 
nents is shown in Fig. 4(b). When the constant current 
pulse is applied (in zero external magnetic field), it in- 
duces initially small-angle in-plane precessions (IPP) of 
the sensing layer around the z-axis. The precessional an- 
gle rapidly increases and spin dynamics turns to the OPP 
regime, where spin precesses mainly around the demagne- 
tization field. The out-of-plane component of the STT, 
Tip, assists in the transition to the OPP-like regimei^. 
Numerical analysis reveals that the transition depends on 
the current amplitude, and spin can precess with positive 
or negative Sx component. Considering constant initial 
spin direction, OPP direction depends mainly on the cur- 
rent density. This appears because the spin phase por- 
trait, and hence the spin trajectory, are modified due to 
the current density. Such a situation is similar to that dis- 
cussed in Fig. 3(c) for the Py/Cu/Py spin valve. Due to 
the sustained energy pumping to the system via the spin- 
transfer, the OPP angle decreases and the spin is finally 
driven into one of the possible static states (SS) close to 
the Ex (SS+), or — e^, (SS_), depending on the sign of 
the Sx component in the OPP regime. The SS± states 
are the static fixed points that result from the interplay 
between the STT and effective magnetic field (mainly 
its demagnetization part). The SS± points are close to 
the maximum magnetic energy. Therefore, if current is 
turned off, the spin position becomes unstable; spin is 
driven due to Gilbert damping through the OPP regime 
with decreasing precessional frequency to the IPP regime. 
In the IPP regime, spin precesses around -t-e^ (— e^) di- 
rection and is finally damped to the P (AP) state. We 
have observed that position of the spin in the SS_ (SS+) 
results in the final P (AP) state, see Fig. 4(b). Thus 
the alternation between P and AP states in the region 
(iv) in the diagram shown in Fig. 4(a), is predominantly 
controlled by the current since the position of the static 
state depends on the current density. 



To elucidate other regions in the diagram [Fig. 4(a)], 
we have to consider shorter pulses. According to the di- 
agram, to have a successful switching event, the pulse 
has to exceed a critical current density and duration. In 
the static limit {tp — > oo) the critical density is about 
/ = — LO/q. For a finite pulse duration, higher densities 
are necessary to drive the spin during the time tp away 
from the P state. When the pulse is shorter than a crit- 
ical one, the relaxation back to the P state takes place 
[region (i)]. To escape the basin of attraction, the time tp 
has to exceed an escape time that depends on the actual 
magnetic energy and the current density. For the higher 
current densities, the spin is driven faster away and a 
shorter escape time is needed. If the pulse ends just be- 
fore the onset of the OPP regime, the spin is then placed 
within the basin of attraction of the AP state. The relax- 
ation via the Gilbert dissipation drives the system to the 
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FIG. 4: (Color online) Current-pulse induced switching in the 
asymmetric Co(8)/Cu(10)/Py(8) spin valve at zero magnetic 
field, (a) Switching time as a function of the pulse dura- 
tion tp and reduced current density ///o, -^o = lO^Acm"^. 
(b) Temporal evolution of s under the 6 ns current pulse of 
the amplitude / = — 3.45 7o. (c) Thermally assisted switching 
probability, Paw, from P to AP, driven by 6ns rectangular 
current pulse, calculated as a function of current density at 
T = 4.2K (solid line) and T = 77K (dotted line). 



AP state [region (ii)]. If the spin is driven further away 
from the P state, the IPP regime switches fast to the 
OPP one. In such a case the final state strongly depends 
on the precession phase at t = tp, i.e., when the current is 
turned off. This gives rise to the P/AP bistability in the 
region (iii). The bistable regime appears up to the pulse 
duration that is not longer than the time necessary for 
spin stabilization in the SS± state. Note, that this anal- 
ysis is valid only for current amplitudes I >3Iq. In case 
of smaller amplitudes, only the steady-state large- angle 
IPP regime has been observed, and apart from the region 
(i) only the region (ii) is present. The periodic rib-like 
structure in the region (ii) arises from the dependence of 
the final state on the precession phase at time t = tp. 

When temperature is nonzero, the final state is af- 
fected by the thermal noise that modifies the overall 
spin switching trajectory. In Fig. 4(c) we show switch- 
ing probability as a function of the current density for 
tp = 6 ns, and T = 4.2K (solid line) and T = 77K (dot- 
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ted line). The switching probabihty oscillates following 
the zero-temperature stripe structure, similarly as for the 
spin valve Py(20)/Cu(10)/Py(8), see Fig. 3(b). However, 
the probability Pg^ oscillates now around the value of 
^sw = 0.5, and for increased current densities approaches 
this value for any temperature. We have found that for 
higher current densities the spin is driven via the OPP- 
like transient regime much closer to the SS state. This 
regime is sensitive to the thermal fluctuations mainly due 
to the component of the thermal field transverse to the 
spin trajectory. When the spin remains in the transient 
regime for a longer time, the impact of the thermal fluc- 
tuations is larger and leads to equilibration of the proba- 
bilities for switching to the P and AP states (Psw 0.5). 



C. Switching in nonstandard spin valves 

The fastest switching process in the asymmetric spin 
valves appears in the region (ii) [see Fig. 4(a)]. This re- 
gion, however, is non-compact and therefore to obtain 
a successful switching one has to set the current pulse 
parameters very precisely. In the region (iii) the bista- 
bility of the final state makes the switching out of con- 
trol. Thus, the most convenient for switching seems to 
be the region (iv). For a proper choice of parameters 
(pulse duration and current amplitude, including also 
the thermal effects), corresponding to the maximum of 
the 

-Pswi see Fig. 4(c), it is possible to obtain controllable 
switching. However, complex spin dynamics, especially 
the ringing which appears after the end of current pulse, 
significantly lengthen the switching time. In practice, 
the longer the switching time is the more sensitive is the 
spin evolution to the external disturbances and tempera- 
ture. Therefore, it is highly desired from the applications 
point of view to shorten the switching time as much as 
possible. Accordingly, we propose here a double-pulse 
switching scheme. The scheme includes two rectangular 
current pulses of certain amplitudes and durations. The 
first pulse of negative current, referred to as destabilizing 
pulse, drives the spin out of its initial position. We note 
that both the coUinear configurations are now unstable. 
The second pulse of positive current, called stabilizing 
pulse, controls the dynamics and drives the spin into the 
final state. Moreover, the stabilizing pulse shortens the 
switching time (suppressing the ringing) via additional 
energy dissipation from the system. 

In Fig. 5(a) we show time evolution of the s^ compo- 
nent due to single current-pulse and double current-pulse. 
Here we consider infinitely long stabilizing pulse that in 
principle has no effect on the spin dynamics for t > 6.5 ns. 
Due to the first current pulse of density I = — 4.0/o, the 
STT drives the spin to the SS state. When the first 
pulse is not followed by the second (stabilizing) one, the 
spin returns back via the OPP and IPP regimes to the 
initial state. In the case of double-pulse, however, the 
stabilizing pulse of / = 2.0 /q drives the spin to the AP 
state. More systematical study reveals that including the 
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FIG. 5: (Color online) Demonstration of double-pulse switch- 
ing scheme making use of second (stabilizing) pulse, (a) Evo- 
lution of Sz spin component. Line (1) corresponds to a 
single tp = 5ns (destabilizing) current pulse of amplitude 
I = — 4.0/0 in zero magnetic field and zero temperature. Line 
(2) shows the evolution of Sz in the case when the first pulse 
is followed by a second (stabilizing) pulse of opposite direction 
and I — 2.0 Iq (the double current pulse is shown by the blue 
dashed line), (b) Thermally assisted switching probability 
Psw from the P to AP states, driven by a 5 ns (destabilizing) 
current pulse (negative) followed by a stabilization current 
pulse of 7 = 2.0 7o and tp —> 00, calculated as a function of 
the reduced current density I/Iq of the destabilizing pulse, 
Jo = lO'^Acm-^ for T = 4.2K (solid line) and T = 77K 
(dotted line). 

stabilizing current pulse of / = 2 /q leads to considerable 
modification of the switching diagram (not shown) . More 
specifically, the region (ii) becomes wider and switch- 
ing times under this current pulses falls down from 4 to 
1 ns. Bistability in the region (iii) becomes reduced, but 
still not completely removed. Finally, in the region (iv) 
the bands related to switching become enlarged, e.g., at 
T = OK, in the range of amplitudes from / ~ — 3.5/o 
to — 4.5J0 one obtains controllable switching for pulses 
tp > 4 ns. Further manipulation of the stabilizing pulse 
amplitude, indeed, enhances overall controllability of the 
switching. 

In addition to the enhanced controllability due to the 
stabilizing pulse, we have observed enhancement of the 
switching probability at finite temperatures. In Fig. 5(b) 
we show the switching probability as a function of the 
current density of 5 ns destabilizing pulse, that is fol- 
lowed by a stabilization current pulse of / = 2.0 Iq and 
tp Qo. For T = 4.2 K and / < -S/q, the proba- 
bility Psvj oscillates with increasing current magnitude. 
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similarly as in the case of a single pulse, see Fig. 4(c). 
The regions of successful switching are then broadened 
and the corresponding amplitude is close to unity. In the 
case of T = 77K, the switching probability in this region 
is roughly constant and approaches Pgw — 0.8. 

IV. SUMMARY AND CONCLUSIONS 

We have studied dynamics of the current-pulse induced 
magnetization switching in both standard and nonstan- 
dard spin valves. The calculated switching diagrams re- 
veal the pulse parameters that are suitable for ultra-fast 
spin switching. We showed that the schemes of optimal 
switching strongly depend on the type of pillar struc- 
ture, and are different for standard and nonstandard spin 
valves. The wavy-like torque in nonstandard spin valves 
introduces stable points in the phase portrait for nega- 
tive currents. These stable points are responsible for the 



P/AP bistability in the pulse switching diagram, which 
hinders obtaining successful switching. Therefore, we 
proposed a switching scheme making use of two current 
pulses, which can efficiently overcome the bistable be- 
havior. Additionally, the application of second, positive, 
current pulse speeds-up the spin dynamics and fasten the 
switching process for several times. We have also shown, 
that the proposed scheme leads to enhanced switching 
probability even at high temperatures. 
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